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Abstract 

We summarize a search for the top quark with the Collider Detector at 
Fermilab (CDF) in a sample of pp collisions at y/s= 1.8 TeV with an integrated 
luminosity of 19.3 pb _1 . We find 12 events consistent with either two W bosons, 
or a W boson and at least one 6 jet. The probability that the measured yield is 
consistent with the background is 0.26%. Though the statistics are too limited 
to establish firmly the existence of the top quark, a natural interpretation of 
the excess is that it is due to ti production. Under this assumption, constrained 
fits to individual events yield a top quark mass of 174 ± I0t\l GeV/c 2 . The tt 
production cross section is measured to be 13.9l 4 \ pb. 

PACS numbers: 14.80.Dq, 13.85.Qk, 13.85.Ni, 02.20.Fh 

The Standard Model has enjoyed outstanding success, yet the top quark, which 
is required as the weak-isospin partner of the bottom quark, has remained unobserved. 
Direct searches at the Fermilab Tevatron have placed a 95% confidence level lower 
limit of M top > 131 GeV/c 2 [TJ. Global fits to precision electroweak measurements 
yield a favored mass of M top = 174^1^ GeV/c 2 §. 
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One expects that, at Tevatron energies, most top quarks are produced in pairs. 
For M top > 85 GeV/c 2 , each top quark decays to a real W boson and a b quark. The 
observed event topology is then determined by the decay mode of the two W bosons. 
About 5% of the time both W bosons decay to eu or \w (the "dilepton mode" ), giving 
two high-Pr leptons with opposite charge, two b jets, and large missing transverse 
energy (JZ! T ) from the undetected neutrinos Q. In another 30% of the cases one W 
boson decays to ev or fiu, and the other to a qq' pair (the "lepton+jets mode"). 
This final state includes a high-P^ charged lepton, JZl T , and jets from the W and the 
two b quarks. The remaining 65% of the final states involve the hadronic decays 
of both W bosons, or the decay of one or both of the W bosons into r leptons. 
These channels have larger backgrounds and are not considered here. This analysis is 
based on a sample of pp collisions at ^fs =1.8 TeV with an integrated luminosity of 
19.3±0.7 pb -1 , collected at the Fermilab Tevatron by the CDF detector|| in 1992-3. 
The details of the analysis are presented in Ref. 

The momenta of charged particles are measured in the central tracking cham- 
ber (CTC), which sits inside a 1.4-T superconducting solenoidal magnet. Outside 
the CTC, electromagnetic and hadronic calorimeters, arranged in a projective tower 
geometry, cover the pseudorapidity region \r)\ < 3.6, allowing reliable measurements 
of the % T . The calorimeters are also used to identify jets and electron candidates. 
Outside the calorimeters, drift chambers in the region \rj\ < 1.0 provide muon identifi- 
cation. A silicon vertex detector (SVX)0, located immediately outside the beampipe, 
provides precise track reconstruction in the plane transverse to the beam, and is used 
to identify secondary vertices that can be produced by b and c quark decays. A 
three-level trigger selects the inclusive electron and muon events used in this analysis. 

In the dilepton search, both leptons are required to have P? > 20 GeV/c and 
to have opposite charge. At least one of the leptons is required to have \r]\ < 1.0 and 
to be isolated 0. In addition, we require $ T > 25 GeV[|?J. To remove background 
from Z production, we reject ee and fifi events with 75 < M« < 105 GeV/c 2 . For 
M top > 120 GeV/c 2 , the two b quarks have significant energy and are detected with 
good efficiency as jets. By requiring two jets with \rj\ < 2.4 and Et > 10 GeVj?|, 
we reduce backgrounds by a factor of four while preserving 84% of the signal for 
M top = 160 GeV/c 2 . To achieve additional rejection against Z — > rr events and 
events with FJ T induced by jet mismeasurement, we require, for FJ T < 50 GeV, that 
the azimuthal angle between the % T and the nearest lepton or jet exceed 20°. No ee 
or /i/i events pass all cuts. Two e/i events survive. 

We use the ISA JET M Monte Carlo program to determine the acceptance and 
the efficiency of the event-selection criteria. The fractional uncertainty in the effi- 
ciency of the two-jet requirement, due mostly to the limited understanding of gluon 
radiation, decreases from 13% for M top =120 GeV/c 2 to 3% for M top =180 GeV/c 2 . 
Other uncertainties in the detection efficiency come from the lepton-identification 
cuts (6%), lepton-isolation cuts (2%), cuts (2%), structure functions (2%), and 
Monte Carlo statistics (3%). The overall acceptance, €dil, for the dilepton search is 
shown in Table [TJ. The number of expected dilepton events from ti production, using 
this acceptance and the theoretical cross section^, is shown in Table 2. 

The dilepton background from WW production is calculated using ISAJET, 
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assuming a total WW cross section of 9.5 pb|K|], and is found to be 0.16±0.06 
events. WW events may contain two jets due to initial-state gluon radiation. The 
treatment of initial-state radiation in the ISAJET calculation is checked using Z+jets 
data, and good agreement is found. The background from Z — > tt is estimated using 
Z — > ee data, where each electron is replaced by a simulated r that decays leptonically. 
This background contributes 0.13±0.04 events. We estimate backgrounds from bb 
and cc using ISAJET to model production processes, and the CLEO Monte Carlo 



program[ll] to model 5-meson decay. The Monte Carlo rates are normalized to a 



sample of e/i data collected with lower trigger thresholds. We estimate 0.10 ± 0.06 
background events from these sources. Backgrounds from hadrons misidentified as 
leptons (0.07±0.05 events) and the Drell-Yan production of lepton pairs (O.lOl^ol) are 
estimated from inclusive-jet and Z data respectively. The total expected background 
is 0.56^o;i3 events, with two candidates observed. 

Events selected for the lepton+jets search are required to have an isolated lepton 
with E T (P T for muons) > 20 GeV and \r/\ < 1.0, and to have Et > 20 GeV0. Events 
containing Z bosons are removed by rejecting events with an ee or /i/i invariant 
mass between 70 and 110 GeV/c 2 . In Table |^ we classify the W candidate events 
according to the multiplicity, Nj et , of jets with E T > 15 GeV and \r}\ < 2.0J7J. The 
dominant background in the lepton+jets search is the direct production of W+jets. 
The ratio of the tt signal to W+jets background can be greatly improved by requiring 
Njet > 3. This requirement has a rejection factor of ^400 against inclusive W 
production while keeping approximately 75% of the tt signal in the lepton+jets mode 
for M top = 160 GeV/c 2 . In the W+ >3-jet sample, we expect 12+2 (6.6+0.7) tt 
events for M top = 160 (180) GeV/c 2 , using the acceptance discussed below and the 
theoretical cross section. We observe 52 events with Nj et > 3. 

The VECBOS Monte Carlo programjl2] can be used to make estimates of direct 
VK+jets production. Table 3 shows the results of a particular calculation which pre- 
dicts 46 events with >3 jets and seven events with >4 jets. The VECBOS predictions 
for > 3 jets have uncertainties of about a factor of two due to the choice of Q 2 scale 
and cannot be used for a reliable absolute background calculation. We have therefore 
developed a technique for estimating backgrounds in the lepton+jets search directly 
from the data. This technique is described below. Other backgrounds (direct bb, Z 
bosons, W pairs, and hadrons misidentified as leptons) contribute 12.2+3.1 events 0. 
Additional background rejection is needed to isolate a possible tt signal. Requiring 
the presence of a b quark, tagged either by a secondary vertex or by a semileptonic 
decay, provides such rejection. 

The lifetime of b hadrons can cause the 6-decay vertex to be measurably dis- 
placed from the pp interaction vertex. When associated with jets with Et > 15 GeV 
and \tj\ < 2.0, SVX tracks with Pt > 2 GeV/c and impact-parameter significance 
\d\/(Jd > 3 are used in a vertex-finding algorithm!]]. Using these tracks, the decay 
length transverse to the beam, L xy , and its uncertainty (typically ai xy ~ 130 /mi) 
are calculated using a three-dimensional fit, with the tracks constrained to originate 
from a common vertex. Jets that have a secondary vertex displaced in the direction 
of the jet, with significance \L xy \/(jL > 3.0, are defined to be "SVX-tagged." 

We use a control sample, enriched in 6-decays, of inclusive electrons (E T > 
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10 GeV) to measure the efficiency for SVX-tagging a semileptonic b jet. We compare 
this efficiency with that predicted by the ISAJET+CLEO bb Monte Carlo and find 
our measured efficiency to be lower than the Monte Carlo prediction by a factor 
of 0.72±0.21. We then determine the efficiency for tagging at least one b jet in 
a ti event with three or more observed jets, et ag , from ti Monte Carlo rescaled by 
the factor determined above. We find e tag = 22 ± 6% independent of top mass for 
M top > 120 GeV/c 2 . The efficiency, esvx, for inclusive ti events to pass the lepton- 
identification, kinematic, and SVX 6-tag requirements is shown in Table 1. The 
number of expected SVX-tagged ti events with Nj et > 3 is shown in Table 2. Six 
SVX-tagged events are observed in the 52-event W+>3-jet sample. 

Rather than rely on Monte Carlo predictions, we estimate directly from our data 
how many tags we would expect in the 52-event sample if it were entirely background. 
We assume that the heavy-quark {b and c) content of jets in VT+jets background 
events is the same as in an inclusive-jet sample Q. This assumption is expected to be 
conservative, since the inclusive-jet sample contains heavy-quark contributions from 
direct production (e.g. gg — ► bb), gluon splitting (where a final-state gluon branches 
into a heavy-quark pair), and flavor excitation (where an initial-state gluon excites a 
heavy quark in the proton or antiproton sea), while heavy quarks in P^+jets back- 
ground events are expected to be produced almost entirely from gluon splitting [TJ 



We apply the tag rates measured in the inclusive-jet sample, parametrized by the 
Et and track multiplicity of each jet, to the jets in the 52 events to yield the total 
expected number of SVX-tagged events from Wbb, Wcc, and fake tags due to track 
mismeasurement. We have tested this technique in a number of control samples and 
use the level of agreement with the number of observed tags to determine the sys- 
tematic uncertainty on the predicted tag rate. The backgrounds from non-W sources 
(direct bb production and hadrons misidentified as leptons) are also determined from 
the data|T]|. The small contributions from Wc, from WW and WZ production, 
and from Z — > rr are estimated from Monte Carlo events. The total estimated 
background to SVX tags in the 52-event sample is 2.3±0.3 events. An alternate 
background estimate, using Monte Carlo calculations of the heavy-quark processes in 
W+jets events and a fake-tag estimate from jet data, predicts a heavy-quark content 
per jet approximately a factor of three lower than in inclusive-jet events and gives an 
overall background estimate a factor of 1.6 lower than the number presented above, 
supporting the conservative nature of our background estimate. 

In the jy+jets sample, the L xy distribution of observed SVX tags is consistent 
with that of heavy-quark jets. The tags in the W events with one and two jets are 
expected to come mainly from sources other than ti decay, and the rate of these tags 
is consistent with the background prediction, with 16 events tagged and 22.1±4.0 
predicted. 

A second technique for tagging b quarks is to search for leptons arising from the 
decays b — > tvX {t = e or //), or b — > c — > £uX. Because these leptons typically have 
lower Pt than leptons from W decays, we refer to them as "soft lepton tags", or SLT. 
We require lepton P T > 2 GeV/c. To keep this analysis statistically independent of 
the dilepton search, leptons that pass the dilepton requirements are not considered 
as SLT candidates. 
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In searching for electrons from b and c decays, each CTC track is extrapolated 
to the calorimeter, and a match is sought to an electromagnetic cluster consistent 
in size, shape, and position with expectations for electron showers. The efficiency of 
the electron selection criteria, excluding isolation cuts, is determined from a sample 
of electron pairs from photon conversions, where the first electron is identified in the 
calorimeter and the second, unbiased, electron is selected using a track-pairing algo- 
rithm. The electron isolation efficiency is determined from ti Monte Carlo events. The 
total efficiencies are (53±3)% and (23±3)% (statistical uncertainties only) for elec- 
trons from b and sequential c decays respectively. To identify muons, track segments in 
the muon chambers are matched to tracks in the CTC. The efficiency for reconstruct- 
ing track segments in the muon chambers is measured to be 96% using J/ip — > 
and Z — ■> decays. This number is combined with the P^-dependent efficiency 

of the track-matching requirements to give an overall efficiency of approximately 85% 
for muons from both b and c decays. 

The acceptance of the SLT analysis for ti events is calculated using the ISAJET 
and CLEO Monte Carlo programs. The efficiency for tagging at least one jet in a 
ti event by detecting an additional lepton with P T > 2 GeV/c is e tag = 16 ± 2%, 
approximately independent of M top . The efficiency, €slt, for inclusive ti events to 
pass the lepton-identification, kinematic, and SLT 6-tag requirements is shown in 
Table 1. The number of expected SLT-tagged ti events is shown in Table 2. We find 
seven SLT-tagged events with Nj et > 3. Three of the seven also have SVX tags. 

The main backgrounds to the SLT search are hadrons misidentified as leptons, 
and Wbb, Wcc production. As in the SVX analysis, we estimate these backgrounds 
from the data by conservatively assuming that the heavy-quark content per jet in 
Vy+jets events is the same as in inclusive-jet events. By studying tracks in such 
events, we measure the probability of misidentifying a hadron as an electron or muon, 
or of tagging a true semileptonic decay. We use these probabilities to predict the 
number of tags in a variety of control samples, and obtain good agreement with 
the number observed. We expect 2.70±0.27 tags in the W+ >3 jet sample from 
these sources. Other sources (direct bb, W / Z pairs, Z — > tt, Wc, and Drell-Yan) 
contribute 0.36±0.09 events, for a total SLT background of 3.1±0.3 events. The 
number of SLT tags in the W+l and VT+2-jet samples, which should have only a 
small contribution from ti, agrees with the background expectation (45 events tagged, 
44±3.4 predicted). Figure [I] shows the combined number of SVX and SLT tags, 
together with the estimated background, as a function of jet multiplicity. 

Each of the analyses presented above shows an excess of events over expected 
backgrounds, as shown in Table [2| The dilepton analysis observes two events with 
a background of 0.56^o;i3. The lepton+jets 6-tag analysis identifies ten events: six 
events with a background of 2.3±0.3 using the SVX tagging algorithm, and seven 
events with a background of 3.1±0.3 using the SLT tagging algorithm, with three of 
these events tagged by both algorithms. For each of these results we calculate the 
probability, V , that the estimated background has fluctuated up to the number of 
candidate events seen or greater. We find Vdil=^%, 'Ps , vx=3.2%, and Vslt=3-8%- 

To calculate the probability V com uned that all three results together are due only 
to an upward fluctuation of the background, we use the observation of 15 "counts" : 
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the two dilepton events, the six SVX tags, and the seven SLT tags. This procedure 
gives extra weight to the double-tagged events, which are approximately six times 
more likely to come from b and c jets than from fakes, and therefore have a sig- 
nificantly smaller background than the single-tagged events. We have checked that 
we understand SVX— SLT correlations by correctly predicting the number of double- 
tagged jets and events in the inclusive-jet sample. We calculate V com uned using a 
Monte Carlo program that generates many samples of 52 background events, with 
fractions of P^+light quark and gluon jets, Wbb, Wcc, and other backgrounds dis- 
tributed according to Poisson statistics with mean values and uncertainties predicted 
by Monte Carlo calculations QXQ . The number of events with heavy-quark jets is scaled 
up to agree with the more conservative background estimate from inclusive-jet data. 
The predicted number of SVX plus SLT-tagged events is obtained by applying the 
measured efficiencies and correlations in the SVX and SLT fake rates. This number is 
combined with a Poisson-distributed number of dilepton background events to deter- 
mine the fraction of experiments with 15 or more counts from background alone. We 
find Vcombined— 0.26%. This corresponds to a 2.8a excess for a Gaussian probability 
function. 

Assuming the excess events to be from tt, we calculate the cross section for tt 
production in pp collisions at y/s=1.8 TeV. The calculation uses the tt acceptance, 
the derived efficiencies for tagging jets in tt events and a revised estimate of the 
background appropriate for a mixture of tt events and background in the 52-event 
lU+jets sample (rather than assuming it to contain all background as above). In 
Tables [I] and ^| we summarize the acceptances, and the theoretical and measured 
cross sections clS db function of M top . 

We have also studied^ distributions of kinematic quantities for the 52 W+>3 jet 
events. If the top quark is very massive the decay jets will typically be more ener- 
getic than jets in W-l-jets background events. One variable with discrimination is 
Et23 = (Et2 + Etz), where Et2 and Etz refer to the E^s of the second- and third- 
most energetic jets in the event. The VECBOS Monte Carlo program predicts that in 



W+jets background events the median of E T23 is 71 GeV, while 93% of HERWIGpi 



ti events (M top = 160 GeV/c 2 ) have E T2 3 > 71 GeV. In the 52-event sample, 39 events 
have E T2 3 > 71 GeV, as do eight of the ten 6-tagged events. This is qualitatively 
consistent with the tt hypothesis; however additional studies in progress are needed to 
reduce systematic uncertainties on the jet energy scale and on the Ex23 distribution 
of the background. 

Assuming that the excess of 6-tagged events is due to tt production, we estimate 
M top using a constrained fit[15| to each tagged event with four jets. Using the 52-event 
W+>3-jet sample, we require a fourth jet with Et > 8 GeV and \rj\ < 2.4. Seven of 
the ten 6-tagged events identified in the lepton+jets analysis pass this requirement. 
These seven events are fitted individually to the hypothesis that three of the jets come 
from one t or t through its decay to Wb, and that the lepton, Et, and the remaining 
jet come from the other t or i decay 0. If the event contains additional jets, only the 
four highest--Er jets are used in the fit. The fit is made for all six jet configurations, 
with the requirement that the tagged jet in the event must be one of the b quarks. 
There are two solutions in each case for the longitudinal momentum of the neutrino, 
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and the one corresponding to the best x 2 is chosen. 

Application of this method to tt Monte Carlo events (M top = 170 GeV/c 2 ) gives 
a distribution with a peak at 168 GeV/c 2 and a rms spread of 23 GeV/c 2 . Fitting 
Monte Carlo W+jets background events to the tt hypothesis yields a mass distribution 
with a broad peak centered at about 140 GeV/c 2 . 

The results of the fits to the seven events are presented in Figure |2]. In this 
sample, 1.4T 2 /" events are expected to come from background jrj. To find the most 
likely top mass from the seven events, we perform a likelihood fit of their mass distri- 
bution to a sum of the expected distributions from PU+jets and a top quark of mass 
M top . The — log(likelihood) distribution from this fit is shown in the inset to Fig- 
ure 2. Systematic uncertainties in this fit arise from the background estimation, the 
effects of gluon radiation on the determination of parton energies, the jet energy scale, 
kinematic bias in the tagging algorithms, and different methods of performing the like- 
lihood fit. Combining these uncertainties yields a top mass of M top = 174 ± lOltil 
GeV/c 2 , where the first uncertainty is statistical and the second is systematic. The 
statistical uncertainty includes the effects of detector resolution and incorrect assign- 
ments of jets to their parent partons. Using the acceptance for this top mass and our 
measured excess over background we find a a {M top = 174 GeV/c 2 ) = 13.9±|;| pb. By 
performing a simple x 2 analysis on the theoretical prediction for the cross section as 
a function of M top , our measured mass, and our measured cross section, we find that 
the three results are compatible at a confidence level of 11% (1.6a). 

We have performed many consistency checks, and have found some features of 
the data that do not support the tt hypothesis. The sample of inclusive Z events 
serves as a control sample for studying the production of a vector boson plus jets, as 
Z bosons are not produced in tt decay. We find two 6-tagged Z+>3 jet events with 
0.64 expected. Both events have four jets and are SVX-tagged. Though the statistics 
are limited, these events could indicate an additional (non-tf) source of vector boson 
plus heavy quark production, not accounted for in our background estimates. Higher- 
statistics checks of the 6-tagging rate in W or Z+l and 2-jet events are consistent 
with expectations. We also find that the measured tt cross section is large enough to 
account for all observed W+4 jet events. The apparent deficit of events from direct 
production of W+4 jets and other backgrounds is a 1.5-2(7 effect. 

Other features do support the tt hypothesis. One of the dilepton candidate 
events is 6-tagged by both the SVX and SLT algorithms, with approximately 0.01 
double-tagged background events (0.13 signal events) expected. This, together with 
the excess of 6-tagged W+jets events, provides evidence for an excess of both Wbb and 
WWbb production, as expected from tt decays. We have performed a kinematic anal- 
ysis of the lepton+jets sample and conclude that it can accommodate the top content 
implied by our measured cross section. Furthermore, a likelihood fit to the top mass 
distributions obtained from the 6-tagged PU+4-jet events prefers the tt+background 
hypothesis over the background-only hypothesis by 2.3 standard deviations. 

In conclusion, the data presented here give evidence for, but do not firmly 
establish, the existence of the top quark. Work is continuing on kinematic analyses 
of the present data, and we hope for an approximate four-fold increase in data from 
the 1994-95 Tevatron collider run. 
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M top 


120 GeV/c 2 


140 GeV/c 2 


160 GeV/c 2 


180 GeV/c 2 


£dil 


0.49 ± .07% 


0.66 ± .07% 


0.78 ± .07% 


0.86 ± .07% 


tsvx 


1.0 ± 0.3% 


1.5 ±0.4% 


1.7 ±0.5% 


1.8 ±0.6% 


tSLT 


0.84 ±0.17% 


1.1 ±0.2% 


1.2 ±0.2% 


1.3 ±0.2% 


( P b) 


38.9+^ u 2 8 


16.9^8 


8 9+1-4 
0-^-0.8 


4 o+U.6 




22.7t^° 9 ° 


16.8^ 


i4.7±|:f 


13.711;? 



Table 1: Summary of top acceptance (including branching ratios) and the theoretical 
cross section 0. The last line gives the it production cross section obtained from this 
measurement. 



Channel: 




Dilepton 


svx 


SLT 


^ expectedi M-top 
^ expectedi -^top 
■^expected-! -^top 
■^expectedi ""-tap 


120 GeV/c 2 
140 GeV/c 2 
160 GeV/c 2 
180 GeV/c 2 


3.7 ±0.6 

2.2 ±0.2 

1.3 ±0.1 
0.68 ±0.06 


7.7 ±2.5 

4.8 ± 1.7 
2.7 ±0.9 
1.4 ±0.4 


6.3 ± 1.3 
3.5 ±0.7 
1.9 ±0.3 
1.1 ±0.2 


Total Background 


0-561^ 


2.3 ±0.3 


3.1 ±0.3 


Observed Events 




2 


6 


7 



Table 2: Number of tt events expected assuming the theoretical cross section, and 
the number of candidate events observed with expected backgrounds. 
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N jet 


Electrons 


Muons 


Total 


Jet 


10,663 


6,264 


16,927 


1 Jet 


1058 


655 


1713 


2 Jets 


191 


90 


281 


3 Jets 


30 


13 


43 


> 4 Jets 


7 


2 


9 



VECBOS (Q 2 =< P T ">) 



157lil 5 7 
267±1? 
39lg 

7+3.2 
' -2.2 



Table 3: Summary of candidate event yields as a function of jet multiplicity. Jets 
have Et > 15 GeV and \r)\ < 2.0. Also shown are the predicted number of W events 
from the VECBOS Monte Carlo program. The uncertainties shown in the VECBOS 
predictions are dominated by the uncertainty in the jet energy scale; the uncertainty 
in the Q 2 -scale is not included. 
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Figure 1: The sum of SVX and SLT tags observed in the J¥+jets data (solid triangles). 
Events tagged by both algorithms are counted twice. The shaded area is the sum of 
the background estimates for SVX and SLT, with its uncertainty. The three-jet and 
four-or-more-jet bins are the tt signal region. 
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Figure 2: Top mass distribution for the data (solid histogram), the W^+jets back- 
ground (dots), and the sum of background + Monte Carlo tt for M top = 175 GeV/c 2 
(dashed). The background distribution has been normalized to the 1.4 background 
events expected in the mass-fit sample. The inset shows the likelihood fit used to 
determine the top mass. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ex/9405005vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ex/9405005vl 
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